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We consider the repulsive Hubbard model on three highly frustrated one-dimensional lattices - 
sawtooth chain and two kagome chains - with completely dispersionless (flat) lowest single-electron 
bands. We construct the complete manifold of exact many-electron ground states at low electron 
fillings and calculate the degeneracy of these states. As a result, we obtain closed-form expressions 
for low-temperature thermodynamic quantities around a particular value of the chemical poten- 
tial fiQ. We discuss specific features of thermodynamic quantities of these ground-state ensembles 
such as residual entropy, an extra low-temperature peak in the specific heat, and the existence of 
ferromagnetism and paramagnetism. We confirm our analytical results by comparison with exact 
diagonalization data for finite systems. 
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I. INTRODUCTION AND MOTIVATION 

The Hubbard model is a particularly simple model 
for strongly interacting electrons in solidsi Neverthe- 
less, rigorous analysis of the model is a difficult task 
and exact statements about its properties are notori- 
ously rare.-i* 2 -' 3 " A number of rigorous /exact results have 
been obtained for the Hubbard model on some specific 

lattices, see, e.g., Refs. miQ!@IIlIII0. In P artic - 
ular, in the context of the origin of ferromagnetism in 
itinerant electron systems^ 3 - different lattices support- 
ing dispersionless (flat) single-electron band were stud- 
ied in some detail^^ ' 14 ' 15 In the last years the the- 
ory of flat-band ferromagnetism has been developed 
further i 16 i 17 ' 18 i 19 i 20 i 21 i 22 Although at first glance one may 
think that the lattices admitting rigorous treatment are 
rather artificial, nowadays new possibilities to design in- 
teracting lattice system with controlled geometry emerge. 
Thus, modern strategies in chemistry open a route to 
synthesize new materials with a desired lattice structure 
and intersite interaction^ Furthermore, recent progress 
in nanotechnology allows the fabrication of quantum dot 
superlattices and quantum wire systems with any type 
of latticed Another rapidly developing field is the con- 
trolled setup of optical lattices for cold atoms . 10 i 24 i 25 

On the other hand, during the past years it has 
been noticed that exact ground states of the quantum 
XXZ Hcisenberg antiferromagnet can be constructed 
at high magnetic fields for a large class of geomet- 
rically frustrated lattices^ These states, called inde- 
pendent or isolated localized magnons, are localized 
on nonoverlapping restricted areas of the lattice and 
they clearly manifest themselves in various peculiari- 
ties of the low-tempcrature strong-field properties of the 



spin systems (macroscopic magnetization jump^ field- 
tuned lattice instability^ 7 - residual entrop y 28 ' 29 ' 30 en- 
hanced magnetocaloric effect ) 28 : 31 order-disorder phase 
transition below the saturation fiel d 29 : 30 etc.). Inter- 
estingly, flat-band ferromagnetism of Hubbard electrons 
exhibits some similarities to the localizcd-magnon ef- 
fect for XXZ Hcisenberg antifcrromagncts on certain 
frustrated lattices! 11 ' 12 : 32 Note, however, that while the 
one-particle description may be identical, the many- 
particle picture is obviously different. While the spin 
model can be viewed as a hard-core bosonic system with 
nearest-neighbor intersite repulsion, the electronic Hub- 
bard model is a two-component fcrmionic system with 
on-site repulsion between different species. Nevertheless, 
it has been found recentl y 11 : 12 ' 32 that several ideas de- 
veloped for the Hcisenberg model can be carried over to 
the Hubbard model. 

In the present paper we consider the repulsive Hubbard 
model on a class of onc-dimcnsional frustrated lattices, 
namely the sawtooth chain and two different kagome 
chains (see Figs.QJi,[Tb, and[Tb). These Hubbard systems 
have highly degenerate ground states for certain electron 
numbers n < n max , n max oc A, where N is the number of 
lattice sites. We will give explicit analytical expressions 
for all these ground states. Their number grows exponen- 
tially with the system size and can be counted by a map- 
ping of the electronic problem onto a one-dimensional 
classical hard-dimcr gas. Moreover, these systems show 
saturated ground-state ferromagnetism for particular val- 
ues of electron number n, i.e., the square of the total 
spin is S 2 = (n/2)(n/2 + 1). Although the sawtooth 
chain on the one hand and the kagome chains on the 
other hand belong to different types of flat-band ferro- 
magnets, the Hubbard model on all three lattices exhibits 
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FIG. 1: (Color online) Three one-dimensional lattices consid- 
ered in this paper: (a) the sawtooth chain, (b) the kagome 
chain I, and (c) the kagome chain II. For the sawtooth Hub- 
bard chain the hopping parameter along the zig-zag path t' 
is v2 times larger than the hopping parameter t > along 
the base line. For the kagome Hubbard chains all the hopping 
parameters t > are identical. Bold (red) lines denote the 
minimal trapping cells for localized electrons. 



an identical thermodynamic behavior at low tempera- 
tures around a certain value of the chemical potential (j® 
if TV — » oo. Indeed, while the sawtooth lattice being a 
one-dimensional version of Tasaki's model is an example 
of the "cell construction" ^ the kagome chain I and the 
kagome chain II belong to Miclke's class of "line graphs" p 
where the kagome chain I is the line graph of the two-leg 
ladde r 33 ' 34 and the kagome chain II is the line graph of a 
decorated two-leg ladder. The number of sites in the unit 
cell is 2, 3, and 5 for the sawtooth chain, kagome chain I, 
and kagome chain II, respectively. More differences can 
be seen in the single-electron energies for these lattices: 
although the tight-binding model of all three lattices ex- 
hibits a dispersionless (flat) lowest-energy band, the next 
dispersive band is separated by a finite gap for the saw- 
tooth chain, but it touches the flat band at one point in 
momentum space in the case of the kagome chains (for 
more details see Sec. IIIip . We mention that all of the 
three lattices were discussed previously in the context of 
vario us pro blems of str ongly c orrelate d sy stems, see, e.g., 
Refs. [35|36|l33l3il3ili3t41 42 43.44 digl 



We consider the standard Hubbard Hamiltonian 
H = fl 0cr + fli nt + (iij^ + n itl ) , 

flint = U y"^»,t n U' ( L1 ) 

i 

The sums run over TV lattice sites i or over the nearest- 
neighbor pairs and periodic boundary conditions 
are imposed, c\ a (ci jCr ) are the usual fermion creation 
(annihilation) operators. For the sawtooth chain the hop- 
ping parameters along the zig-zag path if are y/2 times 
larger than the ones along the base line t > (then 
the lowest single-electron band is completely flat, see 
Sec. IIIip . For the kagome chains all the hopping pa- 
rameters equal t > (then the lowest single-electron 
band is completely flat, see Sec. IIIip . U > is the on- 
site Coulomb repulsion for electrons with different spins. 
With further statistical-mechanics calculation in mind we 
also introduce the term with chemical potential fj,. Note 
that the sign of the term with /i in Eq. (jl.ip is chosen 
to have direct correspondence between the chemical po- 
tential [i for the electronic model and the magnetic field 
h for the respective antiferromagnetic spin mode h 11 ' 12 ! 32 
In what follows we often will set t = 1. 

The remainder of the paper is organized as follows. We 
begin with a brief summary of our main results, Sec. [H] 
In Sec. Mil we discuss the tight-binding model of non- 
interacting electrons [we set U = in Eq. (jl.ip ] for the 
three chains. In Sec. [TV] we construct the complete set of 
exact many-electron ground states of the repulsive Hub- 
bard model on the three lattices for the electron numbers 
n < n max , where n max — TV or TV 4- 1 and TV oc TV, and 
discuss some properties of these states. Moreover, we 
explain the mapping of a certain subset of these states 
onto spatial configurations of classical hard dimers on 
a simple chain. This mapping is crucial to calculate 
the exact degeneracies of the ground states for electron 
numbers n < TV" and to determine the residual entropy 
caused by these states. In Sec. |V] we calculate analyti- 
cally the contribution of the highly degenerate ground- 
state manifold to the grand-canonical partition function 
of the Hubbard model on the respective one-dimensional 
lattices. This contribution dominates at low tempera- 
tures when the chemical potential \i is around a certain 
value /io (fJ>o = 2i for all three lattices). We also cal- 
culate analytically the low-temperature behavior of sev- 
eral thermodynamic quantities such as the average num- 
ber of electrons, the entropy, and the specific heat, and 
we compare these analytical findings with numerical re- 
sults obtained by exact diagonalization for finite lattices. 
Moreover, we use exact diagonalization for finite lattices 
to discuss the influence of small deviations from the ideal 
geometry (leading to a dispersion of the former flat band) 
on the low-temperature thermodynamics. In Sec. I VII we 
use the analytical findings based on the localized-state 
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picture and complementary numerical data from exact 
diagonalization of finite systems for the discussion of the 
ground-state magnetic properties of the considered Hub- 
bard chains. In particular, we discuss the appearance 
of ferromagnetism and paramagnetism. In Sec. I VIII we 
briefly discuss the relation between the electron mod- 
els and corresponding localized-spin models. Finally we 
summarize our findings in Sec. I VILA Some auxiliary cal- 
culations are collected in Appendices. 



II. SUMMARY OF RESULTS 

Compared to previous work on flat-band Heisenberg 
XXZ magnets, the theme of this paper is the new physics 
which arises in itinerant magnets in which (a) the mo- 
bile degrees of freedom are subject to the Pauli principle 
and (b) the Hamiltonian exhibits full SU(2) symmetry; 
in particular, these features set this work apart from pre- 
vious work on XXZ models. 

Previous studies on itinerant flat-band fcrromagncts 
were focused on particular values rif of the electron num- 
ber n for which ground states with saturated ferromag- 
netism exist. By contrast, we characterize the complete 
set of ground states at electron numbers n < ri{, count 
their numbers, and use them to obtain explicitly the low- 
temperature thermodynamic quantities as well as the av- 
erage ground-state magnetic moments. 

Our results cover, on an equal footing, two families of 
fiat-band ferromagnets in one dimension. The first are 
those obtained from Tasaki's cell construction and those 
described by Mielkc's line-graph construction. 

Somewhat unusually for a strongly interacting itiner- 
ant many-body system, for these lattices we provide an 
explicit construction of the full set of exact ground states, 
for a finite range of doping. The construction of the 
ground states is based on a mapping to hard-core dimcrs 
on an appropriate one-dimensional structure. 

This result enables us to obtain the corresponding par- 
tition functions, and hence the low-temperature thermo- 
dynamics, of the magnets in the low-doping regime. 

The most salient consequences concern (i) the entropy, 
(ii) the low- temperature specific heat, (iii) the depen- 
dence of the average number of electrons on the chemical 
potential and (iv) the magnetic properties. 

Regarding (i) the flat band in the one-particle ener- 
gies leads to a huge degeneracy of the many-body ground 
states for a certain range of electron densities resulting 
in a residual entropy in the thermodynamic limit. This 
highly degenerate ground-state manifold has a great im- 
pact on the low-temperature physics in the low-doping 
regime. In particular, an extra low-temperature peak in 
the specific heat appears which is related to an emerging 
low-energy scale separated from the energy scale deter- 
mined by the value of the hopping integral. Moreover, the 
zero-temperature average number of electrons exhibits a 
jump at a certain value of the chemical potential. Re- 
garding (iv) the structure impressed on the many-body 



wave-function by the Pauli principle leads to a degree of 
ferromagnetism in finite systems which varies with elec- 
tron number n; there are specific fillings at which ferro- 
magnetism is saturated [the square of the total spin is 
S 2 = (n/2)(n/2 + 1)] while at others it is only partially 
developed. However, in the thermodynamic limit the re- 
gion of electron density n/N , for which ground-state fer- 
romagnetism exists shrinks to one point (n/N =1/2 for 
the sawtooth chain, n/N = 1/3 for the kagome chain I, 
and n/N = 1/5 for the kagome chain II). For lower elec- 
tron densities the ground state is paramagnetic and the 
low-temperature behavior of the zero-field susceptibility 
follows a Curie law. 

Finally we emphasize here that the localizcd-magnon 
states for the XXZ Heisenberg antifcrromagnct on all 
three lattices can be also mapped onto a one-dimensional 
model of hard dimers. Interestingly, due to the Pauli 
principle, the localized-electron states are less constrain- 
ing than the respective localized-magnon states. As a 
result the manifold of localized states for the electronic 
system is much larger than that for the magnon system. 

III. NON-INTERACTING ELECTRONS. 
TRAPPED ELECTRON STATES 

For zero Coulomb interaction U the diagonalization of 
the Hamiltonian (jl.ip is straightforward. Nevertheless, 
we start with the discussion of this case, since it provides 
some important results which are relevant for the case 
U > 0, too. The sawtooth chain has been studied in 
great detail before (see, e.g., Ref. HH). We will therefore 
just recall the main results for this case and give further 
details only for the two kagome chains. For brevity, we 
may also omit spin indices as irrelevant in this section. 

A. Sawtooth chain 

The sawtooth chain consists of Af = N/2 cells, each 
cell contains two sites, see Fig. [T^,. Hence, there are two 
branches of single-particle energies which read^ 2 - 



£1,2(^0 = t cos k =F y t 2 cos 2 k + 2t' 2 (1 + cos k) + /.t. 

(3.1) 

For t 1 = \[2t > the lowest single-electron band be- 
comes flat, £i(«) = E\ = —2t + fi and one can write down 
creation operators for a corresponding set of eigenstates 
which are localized in a valley with index j££ 

4j = C 2j-1 — ^^ C 2j + C 2j+1 ■ (3-2) 

This set of localized single-electron states is a convenient 
starting point for the construction of the many-electron 
ground states of the Hamiltonian (jl.ip in the subspaces 
with electron numbers n = 2, . . . ,Af. Moreover, this "lo- 
calized" point of view allows a useful simple geometrical 
interpretation, see below. 
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B. Kagome chain I 

Consider next the kagome chain I, see Fig. [TJa. It con- 
sists of A^ = N/3 cells, each cell contains three sites. 
After standard transformations we get the diagonal form 
of the Hamiltonian 



£l(/c) 



EE £ p( K ) a iL a ?v 

p— 1 K 

-2i + n, 

2t cosk + /i, 

2f (1 + cos k) + fi 



(3.3) 



with k = 2irm/Af, m £ Z, -A/"/2 < m < Af/2. The 
lowest energy band is flat, e\(k) = e% = —2t + fx. Note, 
however, that a state with k = ir (it does exist if A/" is 
even) from the dispersive band £2(1*) has also the energy 
Ex = —2t + fx, i.e., the next band touches the lowest flat 
band at k = n. The K-dcpcndent single-electron states 
are given by o^JO), p = 1, 2, 3, with 



Af-i 



E e ~"%> 



l 2,K 



v/2AA(2 + cosk) ^ 

4j _ 4j + l _ 4j+2 + c 3j+3' 

1 V e- iKj iV - J 
^2A? ^ V 3j+1 3j+2 



1 



a 



3 . /v 



2^(2 + cos K ) ^ 



JV-1 



-WSJ 
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3j-2 ^ >-3j-l 



2c 3 ^ + C 3 j +1 + C 3 j +2 



(3.4) 



Owing to the A/"- fold degeneracy of the lowest fiat 
band, one can use alternatively the states Zg - 1 0^ , j = 

0, 1, . . . , Af — 1 instead of the Af states a\ JO). These 
eigenstates are localized states where an electron is 
trapped on a diamond consisting of four sites, 3j, 3j + 1, 
3j + 3, and 3j + 2 (j enumerates these diamond traps 
and varies from to Af — 1). It is easy to check that 

[Hq + fJ-Yli n i>^3j] = £ ^lj- A real-space picture for the 
state with k = n from the dispersive band £2 (ft) is as 
follows: 



4, Jo) = 



N-l 



I2M 



3=0 



3j+l 



-3j+2 



|0). (3.5) 



This eigenstate is not localized on a finite region, but 
the electron is trapped on the two legs [the sites inside 
the strip (i.e., numbers 3j, 3j + 3, etc.) do not appear 



in a); .JO)]. In what follows we call the state a 2 JO) the 
trapped two-leg state. We may also introduce upper- and 



lower-leg states 

4|o) 
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(3.6) 



Obviously, the electron in the eigenstate £j0) (Lj|0)) is 
trapped on the upper (lower) leg. 



C. Kagome chain II 

Finally, we consider the kagome chain II (Fig. [TJ;) 
which consists of Af = N/5 cells, each cell contains five 
sites. Again standard transformations lead to the diago- 
nal form of the tight-binding Hamiltonian 



M E ni 



Ei(k 
e 2 {k 
e 3 (k 
e±(k 

£ 5 (k 



p—1 K 

= -2t + fx, 

= —t\j2 + 2cosk + fx, 

= t — W3 + 2 cos k + fx, 

= tv2 + 2 cos k + fx, 



t + W3 + 2cosk + fx (3.7) 



with k = 2irra/Af, m e Z, -Af /2 < m < Af/2. The state 
with the energy £2(n = 0) (it docs exist for odd and 
even Af) touches the flat band. The K-dependent single- 
electron states are given by <xj, JO), p = 1, 2, 3, 4, 5. For 

the sake of brevity we give a\ JO) for the lowest-energy 
band, only 



f 



M-i 



^(S-cosk) 5r 



-Sj ' \, ■ I 



-5j+2 ~ °5j+3 ' ^5^+4 l '5j+5 



(3.8) 



The corresponding localized states are given by /5JO), 
j = 0, 1, . . . , Af — 1, where an electron is trapped on a 
hexagon consisting of six sites, 5j,5j + l,5j + 3, 5j + 
5, 5j + 4, 5j + 2 (j enumerates these hexagon traps and 
varies from to Af — 1). It is easily verified that [Ho + 
fx rii, = E\l\y A real-space picture for the state 
k = from the dispersive band 62(11) is given by 



4,ol0) = 



1 



AT-l 
2^5 



,t 

'5j+l 



4 



„t 

-5j+3 



-S3+4 



|0). 
(3.9) 

As for the kagome chain I the electron in this eigenstate 
is trapped on the two legs. Again in what follows we 
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call the state a 2 |0) the trapped two-leg state. Again we 
may introduce upper- and lower-leg states 




(3.10) 



Obviously, the electron in the eigenstate L|J0) (L}|0}) is 
trapped on the upper (lower) leg. 

D. Trapped states and destructive interference. 
The geometrical perspective 

It is useful to discuss the appearance of localized elec- 
tron states from a geometrical point of view. As dis- 
cussed above, for all three lattices we can easily single 
out a small area of the lattice which plays the role of a 
"trapping cell" , namely a V- valley for the sawtooth chain, 
a diamond for the kagome chain I, and a hexagon for the 
kagome chain II (see marked regions in Fig. [1]). Solv- 
ing the single-electron problem for the trap one finds the 
lowest-energy eigenfunction oc J2i a « c I|0) with ai = 1, 
a 2 = — \/2, a.3 = 1 (sawtooth chain, the correspond- 
ing energy is —y/2t' < 0), a\ = —ai = 03 = —04 = 1 
(kagome chain I, the corresponding energy is — It < 0), 
or ai = — a2 = 03 = —04 = 05 = — a§ = 1 (kagome 
chain II, the corresponding energy is — 2t < 0). A crucial 
point is that the scheme of the bonds connecting the trap- 
ping cell with the rest should prevent the escape of the 
localized electron from the trap, i.e., the constructed one- 
electron (localized) state should remain an eigenstate of 
the Hamiltonian (|1.1|) on the infinite lattice. It is easy to 
show that a sufficient condition for this is y\ t r ,iaj = 0, 
where the sum runs over all sites i of a trapping cell and 
r is an arbitrary site which does not belong to the trap, 
see also Refs. [26| andfl^. Indeed, the above condition is 
fulfilled if an arbitrary bond belonging to the trap and 
the bonds attached to the two sites of this bond in the 
trap form such a triangle that the electron amplitude 
on all sites outside the trapping cell is zero (destructive 
quantum interference). Note that a similar localization 
mechanism can be caused by a magnetic field for tight- 
binding electrons in two-dimensional structures, where 
the wave packet is bounded in Aharonov-Bohm cages 
due to destructive interference for particular values of 
the magnetic flux4£ 

Interestingly, also the (extended) upper- and lower-leg 
states (|3.6p . (|3.10p fit to this geometrical picture, if we 
interpret each of the two legs as a regular polygon (with 
an even number of sites). Then again two neighboring 
sites of the polygon are surrounded by equilateral trian- 
gles which prevent the electron to escape from the leg. 
Thus, for the kagome chains we have Af localized states 



located on diamonds or hexagons and, in addition, two 
states trapped on the legs (upper-leg state and lower-leg 
state), i.e., in total Af+2 localized states. (Note that for 
the sawtooth chain such additional states do not exist.) 
However, these Af + 2 localized states are not linearly 
independent, since there is one linear relation between 
them, Ll\0)-Ll\0) = \Z2At|0 ), wh ere L\, L\, and A* are 
given in Eq. (|3.6|) or in Eq. (|3.10[) (see also the general 
discussion of linear independence in Appendix [Bj) . As a 
result, there are only + 1 linearly independent local- 
ized single-electron states. These simple arguments are 
in perfect agreement with the more detailed calculations 
presented in Sec. MI Bl and Sec. MI Cl 



IV. TRAPPED ELECTRON GROUND STATES 

FOR U > 

In the previous section, we have found simple highly 
degenerate localized one-particle ground states of nonin- 
teracting spinless electrons which are created by opera- 
tors 4j, 4j> or l h> see E q s - <E3>> 63)' or GS3) I and 

also by a\ „ or a\ operators for the kagome chains, sec 
Eqs. (|3.5p or (|3.9|) ]. It is straightforward to create a set 
of n-electron ground states with 1 < n < Af for the saw- 
tooth chain and 1 < n < Af + 1 for the kagome chains 
by applying n different (i.e., attached to different trap- 
ping cells) operators or ZL [for the kagome chain 

I (II) we may also apply the operator a 2 (a| o)]- The 
energy of these states is nsi, and their degeneracy grows 
exponentially with the system size. 

Now we return to the spinful case of interacting elec- 
trons, i.e., U > in Eq. (jl.ip . Clearly the inclusion of the 
spin does not change the energy for U = but increases 
the degeneracy. Let us denote the degeneracy at U = 
of the ground states of n electrons by g$ (n) . Obviously, 

one has gtf(n) = ( 2 ™™ ax ), n < n max , where n max = Af 
for the sawtooth chain and odd- Af kagome chain I, but 
"max = Af + 1 for even-Af kagome chain I and kagome 
chain II. 

First of all we note that the Hubbard interaction in Eq. 
([l.ip is a positive scmidefinite operator and hence it can 
only increase the eigenvalues of the Hamiltonian (jl.ip . 
On the other hand, among the huge number of localized 
n-electron states being ground states for U = there is 
a considerable fraction of states which do not feel the 
Hubbard interaction term and thus they remain ground 
states with the [/-independent energy ne\ for U > 0. 
However, it is evident that the ground-state degeneracy, 
9M( n )i should decrease as U > is switched on, i.e., 
gj\f{ri) < g$(n). We will consider the cases of the saw- 
tooth chain and of the kagome chains in more detail sep- 
arately. 
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FIG. 2: (Color online) Illustration of a multi-cluster state. 
For each cluster the total spin of the cluster can be flipped 
independently by applying the cluster spin-flip operator S„ i = 

V c f c ■ t 

'j'Gclustcr j,i Ji ' 



A. Sawtooth chain and localized ground states 

Let us recall that for a localized state an electron with 
arbitrary spin a =| , J, is trapped on three contiguous sites 
(V-shaped trapping cell). It is evident that n-electron 
states where the electrons (independently of their spins) 
are located in disconnected V- valleys (i.e., V- valleys with- 
out common sites) are ground states in the n-electron 
subspace having the energy ne\. The explicit expression 
for this type of ground states reads: i| ^ ai - ■ ■ l\j a |0)> 
V|jW.| >2. 

However, these states do not exhaust all ground states 
in the subspace with n electrons. Another type of ground 
states consists of n trapped electrons all with identical 
spin a occupying a cluster of n contiguous V- valleys, e.g., 

Since the interaction term is inactive, this is still an eigen- 
state of the Hamiltonian with eigenvalue ne\. Now 
we take into account the SU(2)-invariance of the Hub- 
bard model fTTjl . i.e., [S~,H] = [S+,H] = [S Z ,H] = 0, 
where 

i 
i 

s * = 5E(4^t-4,^)- ( 4 - 2 ) 

i 

Repeated application of the spin-lowering operator S~ to 
the state (|4.1[) yields 2s + 1 = n + 1 components of the 
spin-(n/2) multiplct 

(S-) m 4. ijT ...4 ( . i+n _ lu |0), m = 0,l,...,n. (4.3) 

Evidently, all these states have the same energy ne\. Us- 
ing the commutation relation 

[S-,ll M ] = ll u , (4.4) 

one can easily write down the explicit expressions for the 
states (|4.3[) in terms of operators l\j only. 

To construct the remaining ground states in the sub- 
space with n electrons with the energy ne± we consider 



now multi-cluster states*^ We consider all possible split- 
tings of n into a sum n = n% + + . . . ; > such that 
clusters of m consecutive occupied V-valleys are sepa- 
rated by at least one empty valley, see Fig. [5] for an ex- 
ample. Starting from a fully spin-polarized multi-cluster 
state we can now independently act with the cluster spin- 
lowering operators, S~ = Ejeduster m c l,l c J : h where 3 
runs over all sites in the cluster rii. This yields many- 
electron ground states which are products of the multi- 
plet components in each cluster. 

In order to compute the ground-state degeneracy 
<?Af( n )j we now need to compute the number Dj^(n) of 
states which are constructed in the manner which we just 
described. This is complicated a bit by the non-trivial 
components of the SU(2)-multiplets consisting of lin- 
ear combinations of products of localized many-electron 
states. However, all coefficients of these linear combina- 
tions are positive and one can choose one state to rep- 
resent the complete linear combination, for example, the 
one where all spins a =| are at the left of each n contigu- 
ously occupied V-valleys and the a =j are at the right 
of the clusters This reduces the counting problem to 
counting the number of configurations Z(n,Af) of three 
states in the trapping cells, namely empty (0) and oc- 
cupied with a =t,|, subject to the constraint that no 
"[-state is allowed to appear as the right neighbor of a 
j-state. This combinatorial problem can be solved di- 
rectly with a 3 x 3 transfer matrix, yielding the canoni- 
cal partition functions Z(n,Af) (see Appendix [A"|) . One 
small additional step then yields Dj^(n): in the sec- 
tor with n = Af and for periodic boundary conditions, 
there are only two allowed configurations according to 
the rules which we just described whereas the SU(2)- 
multiplct has Af + 1 components (Tasaki's ferromagnetic 
ground states). Hence, in the sector with n = J\f and 
for periodic boundary conditions we need to add Af — 1 
configurations by hand. Putting all this together, we ob- 
tain the ground-state degeneracy of the many-electron 
configurations with n < JV electrons localized in Af traps 

gj^ln) = Dtf(n) , 

D*r(n) = Z{n,Af) + (Af -l)5 n ^. (4.5) 

As described above, the last term, (A/ - — 1) S n ^, ensures 
the correct counting in the sector n = Af. Note further- 
more that by writing gj^{n) = Dj^(n) we have implicitly 
assumed that the construction sketched at the beginning 
of this sections yields all ground states. This is indeed 
the case for the sawtooth chain^ see also below. 

An alternative way to compute Z(n,Af) has been de- 
scribed in Ref.Qj]: the combinatorial problem of the three 
states 0, <7 =|, J, in Af traps can be mapped to a hard- 
dimer problem on an auxiliary simple chain with 2Af 
sites. In this mapping, one associates 2 sites to each V- 
valley in order to accommodate either one spin projection 
such that the aforementioned constraint that no f-state 
is allowed to appear as the right neighbor of a |-state and 
the constraint that double-occupancy of a V- valley is for- 
bidden map to hard-dimcr exclusion rules. This auxiliary 
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hard-dimer problem on 27V" sites has been solved with a 
2x2 transfer matrbsii and therefore we may also refer to 
the trapped electron configurations as "hard-dimer con- 
figurations" . Of course, the different ways of computing 
Z(n,Af) are completely equivalent. 

In what follows, we call many-electron ground states, 
which are constructed only from the localized single- 
electron states trapped in a V- valley "hard-dimer states" . 
Their degeneracy Dj^(n) is given by the second line of 
Eq. (|4.5p . The hard-dimer states are the only ground 
states for the sawtooth chain and for the kagomc chain I 
with an odd number of cells TV. For the kagome chain I 
with an even number of cells TV as well as for the kagome 
chain II in addition to the hard-dimer states we have to 
consider also many-electron ground states which involve 
an extended single-electron state with the flat-band en- 
ergy £i - the two-leg state [see Eqs. (|3 . 5|) and (13.91) ]. 



B. Kagome chains and trapped ground states 
involving two-leg states 

Now we consider the ground states in the subspace 
with n = 1,2...,7V+1 electrons for the case when the 
two-leg state comes into play (periodic kagome chain I 
with even number of cells TV or periodic kagome chain 
II). Repeating the arguments elaborated for the sawtooth 
chain, for n = 1, 2, ... ,7V we can construct hard-dimer 
states the number of which is given in Eq. (|4.5p . However, 
the hard-dimer states do not exhaust all many-electron 
ground states for n = 1,2,... ,7V. Indeed, in the sub- 
space with n = 1 electron we have in addition two ground 
states a\ w ^0), u =|, { for the kagome chain I [see Eq. 

(|3.5|) ] or c4o<r|0)' c =TjI f° r tne kagome chain II [sec 
Eq. (I3J1) ]. 

What happens if n > 11 Consider the case n = 2. 
We can construct ground states using the two-leg state 
as follows. Assume both electrons have the same spin 
polarization and the first one is localized within 1 of TV 
(diamond or hexagon) cells whereas the second one is in 
a two-leg state. Obviously this is the ground state with 
the energy 2s\. More states can be generated applying 
S' ± two times. Thus, the degeneracy of the constructed 
ground state is 37V. 

Interestingly, in the subspace with n = 2 electrons 
there is one more possibility to construct an eigenstate 
of the interacting Hamiltonian (jl.ip with the energy 2e\. 
Recall that the electron in the state £T CT |0) [see Eq. (|3.6p 
or Eq. (|3.10p ] is located along the upper leg, whereas the 
electron in the state LjjO) [see Eq. ([376]) or Eq. (j3~T0|) ] is 
located along the lower leg. Therefore, even if two elec- 
trons being in these states have opposite spins, they do 
not feel the Hubbard repulsion (no common sites) and 
the energy of this two-electron state remains 2e\. Thus, 
we find the following (extra) state in the subspace with 
n = 2 electrons for the periodic kagome chains: 



where L) aa , h\ a are defined by Eq. (|3.6p for the periodic 
kagome chain I and by Eq. (|3.10p for the periodic kagome 
chain II. 

Consider next the subspace with n = 3 electrons. 
Again we can assume three electrons to have the same 
spin, two of which arc localized within 2 of TV cells 
whereas the third one is in the two-leg state with the 
flat-band energy. Obviously this is the ground state with 
the energy 3ei. More states can be generated applying 
S three times. Thus, the degeneracy of the constructed 
ground state is 4(^f). Proceeding with such arguments 
for n = 4, ... ,7V + 1 electrons we can easily construct 
the ground states which involve the two-leg state. Obvi- 
ously, it is easy to count their number which is equal to 

(»+!)(£)■ 

We wish to emphasize here that for the kagome chains 
(kagome I with even TV or kagome II chains) we have 
constructed the ground state in the subspace with n = 
TV+1 electrons, which has the energy (TV+l)e_, starting 
from the fully polarized state with TV electrons occupying 
diamond or hexagon trapping cells, respectively, and 1 
electron being in the two-leg state, and then applying 
the S' ± operator TV+1 times. The degeneracy of the 
constructed ground state is TV + 2. 

Summing up, we have arrived at the following formula 
for the degeneracy of the ground states in the subspaces 
with n < TV + 1 electrons for the periodic kagome chain 
I with an even number of cells TV or the periodic kagomc 
chain II: 



9Af{n) = (1 - <5 n ,Af+i) Dtf(n) + (1 - S nfi ) Lj^r(n), 



L N {n) = (n + 1) 



TV 

71-1 



(4.7) 



(extra) = L\ L\ \Q) 



(4.6) 



[n = 0, 1, . . . ,TV for Dj^(n), whereas n = 1, ... ,7V + 1 for 
Ltf(ri)]. 



C. Some properties of the ground states for 

7i ^max 

As we have explained above, the ground states in the 
subspaces with n < n max electrons can be constructed 
either from hard-dimer states only (sawtooth, kagomc 
I with odd TV, n max = TV) or from hard-dimer states 
and the two-leg state (kagome I with even TV, kagome 
II, "max = TV + 1) and the degeneracy of the ground 
state gj^(n) in the subspace with n < n max electrons can 
be determined either according to Eq. (|4.5p (sawtooth, 
kagome I with odd TV) or according to Eq. (|4.7p (kagome I 
with even TV, kagomc II). For the calculation of Z(n,TV), 
see Appendix [A"l 

Fig. [3] illustrates the difference in the ground-state de- 
generacy gjv(n) conditioned by the ground states involv- 
ing two-leg states. We take as examples TV = 7 and 
8. The sawtooth chain never has two- leg state contribu- 
tions, the kagome chain II always has two-leg state con- 
tributions, and for the kagomc chain I the two-leg state 
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FIG. 3: (Color online) Ground-state degeneracies for n — 
1, . . . , n max for the sawtooth chain (Af = 7, empty diamonds; 
N — 8, empty squares), the kagome chain I (TV = 7, empty 
diamonds; Af = 8, filled squares), and the kagome chain II 
(TV = 7, filled diamonds; Af = 8, filled squares). 



£2(71") — £1 = 2t (sawtooth chain) and Ai = (kagome 
chains). Assuming U —* 00 we found for finite systems of 
N = 30 sites A 2 « 0.4044 (sawtooth) but A 2 w 0.0689 
(kagome I) and A 2 ~ 0.0584 (kagome II). These data give 
a hint that the excitation gap for finite kagome chains is 
much smaller than for the sawtooth chain. A finite-size 
extrapolation to N — > 00 for fixed electron density n/N 
suggests a vanishing gap in the thermodynamic limit even 
for the sawtooth chain. Note, however, that in real ma- 
terials with chain structure impurities and other lattice 
imperfections are always present and therefore one deals 
in practice with (an ensemble of) finite chains having a 
gap to the excited states. 

With respect to the contribution of the trapped ground 
states (including the two-leg states) to thermodynamic 
quantities the question arises whether these states for a 
given n < n max are linearly independent. Their linear 
independence can be shown following the lines of Ref. 
IH, for more details see Appendix IBl 



LOW-TEMPERATURE 
THERMODYNAMICS 



contributions appear only for even Af. Consequently, for 
Af = 8, the degeneracies of the two kagome chains are 
identical whereas for Af — 7 the degeneracy of the kagome 
chain I is identical to that of the sawtooth chain. Firstly, 
we observe that the degeneracies increase rapidly with Af 
(note the logarithmic scale of the vertical axis of Fig. [3J . 
Secondly, there is an obvious contribution of the two-leg 
states. 

Next, we discuss some general properties of the con- 
structed many-electron ground states in the subspaces 
with n < n max electrons. As wc have mentioned already, 
these states having the energy net as in the case {7 = 
arc indeed the lowest-energy (ground) states for U > 
since the Hubbard interaction term H lnt = U '^H,'\ n i,i 
in Eq. Ql.ip is a positive scmidefmite operator and can 
only increase energies. 

Another less simple question concerns completeness of 
the constructed ground states. In other words, are the 
constructed ground states the only ground states? Here 
we use numerics for finite systems to check completeness. 
Exact diagonalization data for the sawtooth chain with 
N = 6, 8, 10, 12, 14, 16, the kagome chain I with N = 
9, 12, 15, 18, 21, 24, and the kagome chain II with N = 
15, 20, 25, 30 are in perfect agreement with predictions 
of Sections IIV AI and IIV B[ compare also Fig. [3] These 
numerical findings are supported by general arguments: 
We know that the localized states (including the two-leg 
states) are complete for [7 = 0. Adding a finite U, a 
part of these states will have higher energy, but no new 
ground states will appear, since U can only increase the 
energy. 

We have also estimated the energy gap between the 
ground state and the first excited state A n in the sub- 
space with n electrons. For N — > 00 we have Ai = 



A. Contribution of trapped electron states 

As shown in the preceding section, for all three one- 
dimensional lattices considered the ground states of 
the Hubbard model (jl.ip in the subspaces with n = 
0, 1, . . . , n max electrons have the energy ns\ and the de- 
generacy gj^ij 1 ) [ see Eq s - l)4.5[) . (|4.7|) ]. Now if the chemi- 
cal potential of electrons [i is around fiQ (piQ = 2t for all 
three models) the constructed ground states in the sub- 
spaces with n = 0, 1, . . . , Timax electrons will dominate 
the grand-canonical partition function at low tempera- 
tures due to their huge degeneracy, i.e., 

E(T,n,N) Pa S GS (T,/i,A) 

n max 

= 5Ar(n)exp (-^r) 

n=0 

= m{n) exp 

n=0 L 

^max 

= "£gu(n)z n , (5.1) 

n=0 

where z = exp a; is the activity and x = (/-to — n)/T. 

Consider first the sawtooth chain and the kagome chain 
I with an odd number of trapping cells Af, i.e., all ground 
states correspond to hard-dimer states. Inserting Eq. 
(|4.5p into Eq. (|5.ip we arrive at 

S GS (T, M ,A) = ~ tiap (z,Af) + (Af-l)z Ar , 

Af 

St r ap(^A0 = 5>"Z(n,A0. (5.2) 

n=0 
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Calculating 'E. tYa , p (z,Af) with the help of the transfer- 
matrix method^ (see Appendix [3} we get the following 
equations 

S GS (I>,A0 = # + + g, 



(5.3) 




(Af-l)*z. 



This is an important result, since it allows to calculate 
the contribution of all the ground states described by 
hard dimers to thermodynamic quantities explicitly. 

Consider next the kagome chain I with an even number 
of trapping cells Af and the kagome chain II. Insertion of 
Eq. P~T]) into Eq. (jo~Tj) yields 



S GS (7>, AO = S trap (z,A0 + [Af - l)z 

+ 2^(71+1)1 _ I 



71 = 1 



n- 1 



(5.4) 



Here the first two terms account for the many-electron 
configurations associated to the Af trapping cells, cf. 
Eq. (|5.2p , whereas the third and the fourth term are due 
to the two-leg state. After simple calculations we get the 
following final result 



Sgs(T>,A0 



AT 



i 



(5.5) 



[£± and £3 are defined in Eq. (|5.3[) ]. 

The entropy S{T,n,N) = —dfl(T, fx, N)/dT, 
the grand-canonical specific heat C(T,fi,N) = 
TdS(T,fi,N)/dT, and the average number of electrons 
n(T,fi,N) = dn(T,n,N)/dn follow from Eqs. 
(|5.5|) and the formula for the grand-thcrmodynamical 
potential Cl(T,fj,,N) = -TlnE(T,fj,,N). 

Although explicit formulas for these thermodynamic 
quantities for finite systems are too cumbersome to be 
written down explicitly here, it might be useful to con- 
sider some limiting cases: 

(i) For T = and fi < fj,o we have z — > 00 and con- 
sequently we find for the residual entropy S(fx, N) = 
ln(A/" + 1) and n(fi, N) = Af [hard-dimer states only, cf. 
Eq. (JSU)] or S([x, N) = \n(Af + 2) and n(fx, N) = Af + 1 
[hard-dimer states and two-leg states, cf. Eq. (|5.5[) ]. 

(ii) For T = and /.t > fio we have z — > which leads to 
S((x, N) = and n(fx, N) = 0. 

(iii) At /i = /io we have z = 1 for all T and we find for 
the residual entropy S(m>, N) = ln{[(3 + v / 5)/2] A ^+ [(3 — 
VE)/2] M + Af - 1} [hard-dimer states only, cf. Eq. 1(0]) ] 
or S(fx , N) = ln{[(3 + x/5)/2]^ + [(3 - VE)/2}" + Af + 
2*T+i + j\f2 Af ^ 1 } [hard-dimer states and two-leg states, 
cf. Eq. (|5.5p ]. Moreover, C(T, (J,q,N) = for any temper- 
ature independently of the system size. 

Note, finally, that the conventional (canonical) specific 
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FIG. 4: (Color online) Canonical residual entropy 
S(n,N)/N = \n[Dj^(n)]/Af versus n/Af for system sizes 
Af = 32, 64, 128, 256, 512, 1024. 



heat C(T, n, N) for a fixed number of electrons n < n roax 
is identically zero within the localized-state description. 

In the thermodynamic limit Af — > 00 only the largest 
term, i.e., survives in the expressions (|5.3p and t|5.5|) 
for the partition function^ S G s(T, A* ; N) and we have 



E GS (T,fx,N) 




(5.6) 



which holds for all three considered one-dimensional lat- 
tices. The only differences consist in the relation be- 
tween the number of cells Af and the lattice size N 
(Af = N/2, Af = N/3, and Af = N/5 for the saw- 
tooth chain, the kagome chain I, and the kagome chain 
II, respectively). Explicit formulas for these thermody- 
namic quantities can be obtained easily from the sim- 
ple expression (|5.6p for S G g(T, /x, N). For instance, the 
residual entropy for Af —> 00 is given by S((xq, N)/Af = 
ln[(3 + V5)/2] = 0.962 42 . . . , i.e., the ground-state de- 
generacy grows asymptotically with the system size ac- 
cording to tp 2J ^ , where interestingly the golden mean 
ip = (1 + v / 5)/2 enters the expression. 

Using the transfer-matrix representation we can also 
calculate numerically the (canonical) residual entropy 
S(n,Af)/Af = ln{Dtf(n)}/Af related to the hard-dimer 
states for certain fixed values of the electron number n 
for Af up to 1024. The results are shown in Fig- [4j From 
Fig. 2] it is obvious that the largest degeneracy is found 
for n around Af/2. The extrapolation to Af — > 00 yields 
limv-»oo ln[Djv-(n = Af/2)]/Af « 0.955, which is very 
close to the grand-canonical residual entropy at fi = fxo, 
see above. Indeed, the maximum of limj^— >oc S(n,Af)/Af 
should reproduce the grand-canonical residual entropy, 
albeit with slow finite-size convergence. 
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FIG. 5: (Color online) Average number of electrons in the 
ground state n(0, /x, N) /N versus chemical potential /x//xo for 
the sawtooth chain [N = 12 (solid), N = 20 (long-dashed)], 
the kagome chain I [N = 18 (short-dashed)], the kagome 
chain II [N = 20 (dotted)] with U -y oo. The result for 
n(0, /x, N) /Af which follows from Eqs. (|5.3|l and (|5.5[) is given 
by 0(1 - /x//xo) and [(Af + 1)/Af|0(l - /x//x ), respectively. 



B. Comparison with exact diagonalization 

The above predictions are expected to be valid at low 
temperatures and for /x around /xo. In order to verify 
these expectations and examine the region of validity 
more precisely, we have performed complementary exact 
diagonalization of the full Hubbard Hamiltonian (|l.lj) . 
First, we perform a symmetry reduction of the problem. 
In particular, we use separate conservation of the num- 
ber of electrons with a given spin projection n-\ and nj. as 
well as translational invariance. Ground states can then 
be obtained with the help of the Lanczos recursiom 51 i 52 
In order to obtain thermodynamic properties, we per- 
form a full diagonalization in each symmetry subspace 
using a library routine. The problem simplifies a bit in 
the limit U — > oo where states with doubly occupied sites 
can be eliminated from the Hilbert space. Nevertheless, 
a full diagonalization of the complete problem can only 
be performed for small N . Low-temperature approxima- 
tions for somewhat bigger N can be obtained by omitting 
certain sectors of n. 

Fig.[5]shows exact diagonalization data for the average 
number of electrons per cell n(T, fx, N) /Af at T = as a 
function of the chemical potential /jl/hq for the sawtooth 
chain (N = 12, 20), the kagome chain I (N = 18), and 
the kagome chain II (N = 20). 

In Figs. [6] and [7] we show exact diagonalization data 
for the temperature dependence of the entropy per 
cell S(T,fi,N)/M at /x = 0.98/i , Mo, 1.02/xo [pan- 
els (a)] and the temperature dependence of the grand- 
canonical specific heat per cell C(T, /x, N)/Af at /i = 
0.98/io, jUoi 1-02/xo [panels (b)] for the sawtooth chain 



(N = 12) and the kagome chain I (N = 9) (Fig. EJ) and 
for the kagome chain I (N — 12) and the kagome chain II 
(N = 15) (Fig. [7j) with U — ► oo. We also report analyti- 
cal predictions for finite Af [obtained from Eqs. (|5.3[) and 
(|5.5p ] as well for infinite Af [obtained from Eq. (|5.6p ]. 

From Fig. [5]we see that in the ground state the average 
number of electrons per cell n(T, /i, N) /Af as a function of 
n/fj-o exhibits a jump at jU/yUo = 1. This jump is related 
to the fact that all ground states with < n < rt max 
are degenerate at /i = /io. For — > oo the jump is of 
height unity. For fi/fio < 1 a plateau appears, see Fig. 
[5] The plateau width (i.e., the charge gap) increases as 
U increases and rapidly approaches a saturation value 
for all three finite lattices (for the finite sawtooth chains 
with N = 12, 16, 20 such data were reported in Fig. lb 
of Ref. [lit). What happens with the plateau width as 
increases? For the sawtooth chain the plateau width is 
almost independent of A^ (compare the results for A^ = 
12, 16, 20 in Fig. lb of Ref. [n|). For U -> oo we can 
also find the exact ground state for n = n roax + 1, see 
Sec. I VI CI which allows to determine the size- independent 
plateau width A/x = 2i. By contrast, for the kagome 
chains the plateau disappears as A^ — > oo. 

From Figs. [6] and [7] (temperature dependences of en- 
tropy and specific heat) we see that analytical predictions 
for finite Af as they follow from Eqs. (|5 .3|) and (|5.5p per- 
fectly reproduce the exact diagonalization data at low 
temperatures. Note, however, that for finite N the devi- 
ation from the hard-dimer description is noticeable, al- 
though in the thermodynamic limit A^ — > oo Eqs. (|5.3p 
and (|5.5p imply the one-dimensional hard-dimer behav- 
ior. Finite-size effects are clearly seen in the panels corre- 
sponding to /x = 0.98/xo (Figs.[B]and[7]) and to /x = and 
fi = 1.02yUo (Fig- EI)- The most prominent features seen 
in these plots are a finite value of the entropy at very low 
temperatures for /x = (1q and an extra low-temperature 
maximum in the specific heat C(T, fi, N) for /x 7^ /io- 
The value of the residual entropy for finite systems as it 
follows from analytical predictions based on Eqs. (|5.3p . 
(|5.5p and exact diagonalization are in perfect agreement. 
A high-temperature maximum of C(T, /x, N)/Af around 
T w Mo is common for any system with a finite band- 
width, whereas a low-temperature peak of C(T, /x, N)/Af 
around T w 0.01/xo for |/x — (1q\ = 0.02/xo emerges due 
to the manifold of localized electron states. Indeed, the 
analytical predictions which follow from Eqs. (|5.3p . (|5.5p 
for finite Af are indistinguishable from the U — > 00 ex- 
act diagonalization data in the low-temperature peak re- 
gion. We should, however, mention that for A" — > 00 
the excitation spectrum of the full Hubbard model in the 
many-electron sectors is most likely gapless, in particular 
for finite values of U (compare Sec. II V Cl and Ref. l37lb). 
Such gapless excitations could give rise to quantitative 
corrections at all temperatures for fi < /xq. Indeed such 
quantitative deviations are visible for instance on finite- 
size sawtooth chains with /x = 0.98 /xo and U = 4 ti 11 ^ 7 

To demonstrate the effect of finite U on temperature 
dependences we consider as an example the kagome chain 
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FIG. 6: (Color online) Entropy S(T, fi, N)/Af and grand- 
canonical specific heat C(T, /i, N) /N for the sawtooth chain 
(N = 12, triangles) and the kagome chains I (N — 9, di- 
amonds) with t = 1, U — ► oo. (a) S(T, /j,, N)/Af versus 
temperature T/pn at /i = 0.98/io, /Uo, 1.02/io (from top 
to bottom), (b) C(T, fx, N)/Af versus temperature T/fio at 
/i = 0.98/io, Mo, 1.02/io (from top to bottom). We also show 
the hard-dimer result for N — ► oo as it follows from (|5.6|l (thin 
solid lines) as well as the results which follow from Eq. (|5.3[) 
for M = 3 (dotted lines) and N = & (dashed lines). Note that 
for these systems no additional leg states exist. Note further 
that for fj, = fig and /i — 1.02/io the hard-dimer results for 
J\f = 3, 6, and oo are indistinguishable. 



FIG. 7: (Color online) Entropy S(T,fi,N)/Af and grand- 
canonical specific heat C(T, /x, N) /M for the kagome chain 
I (N = 12, diamonds) and the kagome chain II (N = 15, 
the sectors with up to 7 electrons were taken into account, 
pentagons) with t = 1, U — > oo. (a) S(T, [i, N) /AT ver- 
sus temperature T//io at /i = 0.98/io, /Uo, 1.02/io (from top 
to bottom), (b) C(T, n, N) /N versus temperature T/fio at 
/i = 0.98/io, no, 1.02/io (from top to bottom). We also show 
the hard-dimer result for N — > oo as it follows from (|5.6|l (thin 
solid lines) as well as the finite-size results for M = 3 (dash- 
dotted lines) and M = 4 (dotted lines) which follow from Eq. 
(|5.5p and include the contribution of the leg states. 
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FIG. 8: (Color online) C(T, fj,, N)/N versus T/> for the 
kagome chain I with N = 12 sites and t = 1 for different 
values of U [U = 1 (empty circles) , U = 10 (filled diamonds) , 
(7 — > oo (empty diamonds)]. The results for finite U were 
obtained taking into account the sectors with up to 8 elec- 
trons. We also show the result which follows from Eq. (|5.5p 
for N = 4 (lines) . 



I with N = 12 sites and show C(T, ^, N)/Af versus T//x 
for U = 1, 10, oo (see Fig. [8]). For finite U the calcu- 
lation of thermodynamic quantities becomes very time 
consuming and therefore we report the contribution of 
the subspaces with up to a certain number of electrons 
which is less than 2N thus restricting these data to not 
too high temperatures. As can be seen from these data, 
the features at sufficiently low temperature do not de- 
pend on the value of U > 0. Note, however, that at 
[7 = 1 there are some visible corrections down to tem- 
peratures T/> = e>(icr 2 ). 

From the experimental point of view it is important 
to discuss the stability of the features determined by lo- 
calized electron states with respect to small deviations 
from the ideal lattice geometry Then the conditions for 
the existence of flat bands are violated, and, as a result, 
the exact degeneracy of the ground states in the sub- 
spaces with n = 1, . . . , n max electrons is lifted. However, 
we are still faced with a set of a large number of low- 
lying energy levels which may dominate low-temperature 
thermodynamics as fi is around fiQ. Deviations from ideal 
geometry for the sawtooth chain (i.e., when t 1 slightly de- 
parts from the value \[2t) has been discussed already in 
Refs.HUb and [37lb . Below we consider the kagome chain I 
assuming the hoppings along legs to be slightly different 
from the hoppings along diamonds. Specifically, we put 
t t i = 1 along the diamonds and = 1.01 along the two 



N=24,t L =l, t d =\, n=A 




FIG. 9: (Color online) Canonical specific heat C(T,n,N)/n 
versus T for ideal (td = th = 1) and distorted (td = 1 and 
tL = 1.01) kagome chains I with n — N/6 electrons; U — > 00. 



legs. In Fig. [S] we report the results for the temperature 
dependence of the canonical specific heat C(T,n,N)/n 
for n = Af/2 = N/6 electrons. For the ideal lattice 
C(T,n, N) = if n < n max in the low-temperature 
regime (Fig. [5]). A deviation from exact degeneracy 
{td = 1 7^ tL = 1.01) produces a low-temperature peak 
(Fig. [9]). This peak indicates a separation of two en- 
ergy scales, one is related to the manifold of low- lying 
trapped states, and the other one to the ordinary ex- 
tended states. The position of this extra peak depends 
on td — tL, for instance for td = 1 7^ = 1.1 and N = 24 
the peak is at T = 0.045. Although there are finite-size 
effects in the height and the position of the peak its exis- 
tence is not questioned. We notice that the temperature 
dependence of the specific heat C(T,n, N)/n is an ex- 
perimentally accessible quantity and its well-pronounced 
low-temperature features conditioned by localized elec- 
tron states for small deviations from ideal lattice geom- 
etry may increase chances to observe localized electron 
state effects. 



VI. MAGNETIC PROPERTIES 

The study of ground-state magnetic properties of the 
Hubbard model (jl.lj) on the considered lattices is of great 
interest and has been discussed since the early 1990s. 
Thus, the sawtooth chain belongs to the one-dimensional 
version of Tasaki's lattice. H. Tasaki proved that the 
ground state of the sawtooth-chain Hubbard model is 
ferromagnetic and unique [up to the trivial SU(2) degen- 
eracy] if the number of electrons n = N/2 (saturated 
ground-state ferromagnetism) , see Ref. 0a. Later on, nu- 
merical studies of sawtooth chains of up to N — 12 sites 
by Y. Watanabe and S. Miyashita^ revealed ground- 
state ferromagnetism (saturated and nonsaturated) for 
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other values of n. Moreover, it was shown in Ref. [ll| that 
within the localized electron picture the model exhibits 
full polarization in the ground state for n = N/2 — 1 elec- 
trons (only one-cluster states constitute the set of ground 
states) and 60% of the full polarization in the ground 
state for n = N/2 — 2 electrons if N — > oo. 

The two considered kagomc chains are line graphs. 
This has been pointed out explicitly for the kagome chain 
1 41 where this connection has been known for a while^ 
albeit using a different terminology. According to a gen- 
eral theory elaborated by A. Mielke the repulsive Hub- 
bard model on these lattices should have ferromagnetic 
ground states for the number of electrons n < N / 3 [peri- 
odic odd- Af kagome chain I (nonbipartite parent graph)] , 
n < N/3 + 1 [periodic even- Af kagome chain I (bipartite 
parent graph)], or n < N/5 + 1 [periodic kagome chain II 
(bipartite parent graph)] (Theorem of Ref. 0a or Theo- 
rem 1 of Ref. 03). Moreover, according to A. Mielke the 
ferromagnetic ground state is unique apart from degener- 
acy due to SU(2) invariance (i.e., saturated ground-state 
ferromagnetism) , if n = N/3 (periodic odd- A/" kagomc 
chain I), n = N/3 + 1 (periodic even- Af kagomc chain I), 
or n = N/5 + 1 (periodic kagome chain II) (Theorem 2 
of Ref. [lb). 

To reveal the ferromagnetic ground states we consider 
the operator 



Af 2 



| (S+ S~ + S~ S+) + S z 

JT 2 



(6.1) 



with the operators S a defined in (|4.2|) . The average at 
T = 0, (S 2 ) n /Af 2 is given by the equal- weight average 
over all degenerate ground states for the given number of 
electrons n. It satisfies 



< 



Af 2 



< 



i) 



Af 2 
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2' 



(6.2) 



If (S 2 ) n /Af 2 achieves its maximal value the ground state 
in the subspace with n electrons is the saturated fer- 
romagnetic ground state. If (S 2 ) n /Af 2 has a nonzero 
value which is less than the maximal value the ground 
state in the subspace with n electrons contains ferromag- 
netic ones and is a nonsaturated ferromagnetic ground 
state. To examine ground-state magnetism - in particu- 
lar the existence of ferromagnetism - for thermodynami- 
cally large systems one has to consider the limit A^ — > oo , 
n — > oo preserving n/Af = const. 

Below we use the constructed many-electron ground 
states for n = l,...,n max to discuss systematically 
ground-state ferromagnetism at electron densities n < 
n max . We complete our analytical arguments by numer- 
ics for finite systems. Our results are consistent with 
general theorems of Tasaki and Mielke for n = n roax and 
go beyond considering n < n max . Moreover, we also re- 
port numerics for higher electron densities n > n max . 

We start with a brief overview of our findings: (i) All 
ground states are fully polarized (saturated ferromag- 
netism) for n = Af and n — Af — 1 (all chains) and, 



in addition, for n = Af + 1 for the kagome chain II and 
the even- Af kagome chain I. (ii) For smaller n < Af — 1 
the ground-state manifold contains fully and partially po- 
larized as well as paramagnetic states. Thus we have 
{S 2 ) n < <S , max (5' max + 1). The magnetic polarization 
(S 2 ) n < Smax^max + 1) decays monotonically with in- 
creasing Af — n, n < Af, sec Figs. [TO] and [TTJ For large Af 
the decay becomes very rapid, (iii) While for finite sys- 
tems there is a finite region of electron density n/Af where 
ground-state ferromagnetism exists, this region shrinks 
to one parameter point n/Af = 1 for Af — > oo. (iv) For 
n/Af < 1 the system shows Curie-like behavior with a 
uniform zero-field magnetic susceptibility x ^ T^ 1 . 

In what follows, we first discuss separately the case 
of the systems with hard-dimcr ground states only (see 
Sec. I VI A)) and the case of the systems with ground states 
which also involve two- leg states (see Sec. IVI B[) . In 
Sec. IVI Cl we report exact diagonalization data for finite 
systems at higher electron densities n > n max . We com- 
plete our discussion considering the low-temperature be- 
havior of the uniform zero-field magnetic susceptibility 
for n < n max in Sec. IVI Dl 



A. Hard-dimer ground states and ground-state 
magnetism 

To calculate (S 2 ) n we may further elaborate the 3x3 
transfer-matrix technique. First, we use SU(2) invari- 
ance of the Hubbard model (jl.ip to write (S 2 ) n = 
3{S z2 ) n . Second, we notice that the average over all 
degenerate ground states for a given number of elec- 



trons n is (S z ), 



where S z is the 



z-componcnt of the spin operator of the trap j. The op- 
erator Sj acting on the hard-dimer ground states yields 
(empty trap), 1/2 (occupied trap with spin- up electron), 
— 1/2 (occupied trap with spin-down electron) . Next with 
the help of the 3x3 transfer matrix we find the grand- 
canonical zz correlation function (SqS z ) z (the subscript 
z denotes the activity) yielding the required canonical zz 
correlation function (SQS z ) n (see Appendix A). As a re- 
sult we obtain (S 2 ) n for n = 1, . . . , Af— 1. We recall that 
(S 2 )x = (Af/2)(Af /2 + 1). An alternative computation 
based on a hard-dimer mapping (sec Appendix A) yields 



the same results for (S 



■ ,Af. 



Using the transfer-matrix method we have calculated 
{S 2 ) n /Af 2 , n = l,...,Af for systems of up to Af = 256, 
see Fig. [TUJ Wc obviously have (S 2 ) n = S max (S max + 1) 
for n = Af and Af — 1 independently of the system size 
as predicted above. From the data shown in Fig. [TO] 
we see clearly that (S 2 ) n /Af 2 becomes smaller for any 
fixed n/Af, < n/Af < 1, as Af increases (compare 
the curves for Af = 8 and Af = 256). More precisely, 
the curves shown in Fig. [TO] suggest {S 2 } n /Af 2 — » 
for Af — ► oo at n/Af = const. Note that for satu- 
rated ferromagnetic ground states (S 2 ) n /Af 2 would equal 
(n/Af) 2 / A in the thermodynamic limit Af — > oo. In 
Fig. [TOJ one observes increasing differences between the 
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FIG. 10: (Color online) Average ground-state magnetic mo- 
ment of the sawtooth-Hubbard chain: (S 2 )„/Af 2 versus n/Af 
for Af = 8, 16, 32, 64, 128, 256. For saturated ferromagnetic 
ground states {S 2 )„/Af 2 would equal (n/A/") 2 /4 (thin line) in 



A* 




the thermodynamic limit Af — > oo. 
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FIG. 11: (Color online) Average ground-state magnetic mo- 
ment of the kagome-Hubbard chains: (S 2 ) n /Af 2 versus n/Af 
for AT — 8, 16, 32, 64, 128, 256. For saturated ferromagnetic 
ground states (S 2 )„/Af 2 would equal (n/Af) 2 /^ (thin line) in 
the thermodynamic limit Af — » oo. 



TABLE I: Ground-state {S 2 ) n for n = 1, . 
the sawtooth-Hubbard chain. 
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curve (n/Af) 2 /4 and the data for (S 2 ) n /Af 2 which ap- 
proaches zero for n/Af < 1. We have performed a finite- 
size extrapolation of (S 2 ) n /N 2 for several fixed electron 
densities n/Af = 1/2, 3/4, 7/8 to estimate the cor- 
responding values for Af — > oo, and we find, indeed, 
limvv'-+oo('S' 2 )n/A/' 2 — for those values of n/Af. 



In a next step we will give analytical predictions for 
(S 2 ) n /Af 2 . We first collect the values for (S 2 ) n for sys- 
tem sizes up to Af — 8 in Table fl] Note that the exact- 
diagonalization data for (S 2 ) n for the sawtooth chain 
confirm the values given in Table |TJ By inspecting the 



values of (S 



(S 2 ) 



reported in Table U we may notice that 
3(A/"-3) j\f- 



(S 2 ' 



3 — 
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16(^-6) 

This leads to a guess for the thermodynamic limit: 



(6.3) 



(6.4) 



This result is in accordance with the numerical results 
presented in Fig. [TD1 and gives again evidence for a para- 
magnetic ground state in the infinitely large system. 

Summarizing the above analysis we conclude that there 
is no finite range of ground-state ferromagnetism for elec- 
tron densities n/Af < 1 as Af — » oo. 



B. Ground states involving two- leg states and 
ground-state magnetism 

We turn now to the periodic kagome chain I with an 
even number of cells Af and the kagome chain II for the 
number of electrons n < Af + 1. As was explained in 
Sec. IIV1 the states involving two-leg states increase the 
ground-state degeneracy for n = 1, . . . ,Af by the number 
Lj\f(n) [see Eq. (|4.7[) ]. All these additional states [except 
the one state for n = 2, sec Eq. (|4.6p ] are fully polarized, 
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i.e., S 2 — (n/2)(n/2 + 1). Therefore, for the average 
value we have (S 2 ) n > (S 2 ) n \n, where {S 2 ) n \£> is the 
corresponding value considering hard-dimer states only 
(i.e., it is that value considered in Sec. IVI A[ see, e.g., 
Table [J). Recalling that D_^(n) is the number of hard- 
dimer ground states one finds 



(S 2 ) n = 



D M {n) 



In D 



Dx(n) + Lx(n) 



n /n 
2 V2 



(6.5) 



This formula is valid for n = 1 , . . . , Af — 1 . We re- 
call that {S 2 )m = (Af/2)(Af/2 + 1) and (S 2 )m+i = 
[(Af + 1)/2}[(JV + l)/2 + 1]. Based on Eq. dB3]) we 
again can calculate (S 2 ) n for the kagome chains of up 
to Af = 256 cells, see Fig. [TT] Noticeable deviations from 
(S 2 ) n \D appear only for small Af. Hence the data shown 
in Fig. [TT] give clear evidence that the contribution of the 
two-leg states becomes irrelevant for large systems. This 
conclusion is supported by further inspection of Eq. (|6.5|) . 
Using Eq. (|6.4|) to replace (S i2 )„|_d in the limit AT — > oo 
we can write Eq. (|6.5p for large Af as 



N 



lim (S 2 ) n = R n 



1 



Rn = lim 



„ n n in 
D N (n) 



(6.6) 



According to Eqs. (|A3|) . (|A4j) . ([53]) . and ([53]) the quan- 
tity Dj^f(n) is proportional to the nth derivative of 
+^ + ^3 / " with respect to z at z = 0, whereas L_\f(n) 
is proportional to the nth derivative of t;-^ +£,^ with 
respect to z at z = 0. Evaluating derivatives at z — we 
find that for Af — > 00 the ratio L_\f(n) / D_\f(n) cx 1 /A/" and 
lim J v_too Lj\f(n)/Dj\f(n) = 0. Evidently, i?„ = 1 in Eq. 
(|6.6p implies that Eq. (|6.4p is valid again, i.e., the ground 
state is paramagnetic for the infinitely large system. 



C. Results for n > n max 

The flat-band ferromagncts discussed above may ex- 
hibit ferromagnetic ground-state ordering even for elec- 
tron numbers n > n max . However, this ground-state 
ferromagnetism occurs only for sufficiently large U and 
therefore is definitely different from the true flat-band 
ferromagnetism which emerges for any arbitrary small 
U. 

The occurrence of a saturated ferromagnetic ground 
state can be well understood for the sawtooth chain with 
n = n max + 1, ima = Af. For U = the sawtooth 
chain has two single-electron bands separated by the en- 
ergy gap Ai = £2(71") — E\ = 2t, see Eq. (|3.1[) . If U is 
small (in comparison with, e.g., Ai) the ground state 
in the subspace with n = Af + 1 electrons is a compli- 
cated many-body state. However, if U is sufficiently large 
U > U c (Af+ 1) (and in particular in the limit U — > 00) it 



might be energetically favorable to avoid Hubbard repul- 
sion. This can be realized by occupying all trapping cells 
(V-valleys) with A/", say, spin-up electrons and putting 
the one remaining electron also with a =f into the next 
(dispersive) band. Indeed, it is easy to show that for the 
periodic even-A/" sawtooth chain such a state 



IW+i) 



Lt 2,7r,T < '0,T 2,T 



it 



>T |0) 



is a true eigenstate with the energy Afei + £2(71") 
trapping-cell operators l\ - ^ are defined in Eq, 



(6.7) 

The 
and 

a 2 it t crea tes an electron in the dispersive band with n = 
it and a =J. Other states belonging to a spin- [(Af + 
l)/2] SU(2) multiplet can be obtained by applying S~ 
to the state (|6.7p . This kind of saturated ground-state 
ferromagnetism for the sawtooth chain was first found 
numerically by Y. Watanabe and S. Miyashita— (see also 
Refs.0112). 

We may expect such fully polarized ferromagnetic 
ground states for sufficiently large U for further electron 
numbers n max < n < N, where n = N — 1 is that electron 
number where the well-known Nagaoka theorem^ holds. 
Using Lanczos exact diagonalization of finite systems we 
have investigated this question for the sawtooth and the 
kagome chains I and II. We list our numerical findings 
for U — > 00 in Table [Til Indeed, fully polarized ferromag- 
netic ground states exist for various electron numbers n 
in the range n max < n < N . For other values not listed 
in Table [Til the ground state is either partially polarized, 



< (S 2 ) n < 



1), or it is a singlet, i.e. 



i.e. 

(S 2 ) n = 0, with spiral structure, see also Ref. 1 191. How- 
ever, a detailed discussion of this issue goes beyond the 
scope of the present paper. 



D. Low-temperature behavior of the magnetic 
susceptibility 

We complete our discussion of the magnetic proper- 
ties with a brief consideration of the low-temperature 
behavior of the uniform zero-field magnetic susceptibility 
X- Using the standard arguments for deriving the uni- 
form zero-field Langevin susceptibility we may write the 
trapped-state contribution to \ as 



x(T,n,N) 



3T 



(6.8) 



with {S 2 ) n calculated in Sees. IVI AI and IVI Bl Thus we 
may expect a Curie-like behavior of the susceptibility of 
the considered Hubbard chains at low temperatures in 
case of paramagnetic ground states. In the thermody- 
namic limit we have lim./v'^oo(<S' 2 )n = 3n/4 in the para- 
magnetic region (n/Af < 1), see Eq. (|6.4p . Therefore 
the Curie constant is n/4, which corresponds to a sys- 
tem of n independent spins 1/2. Note, however, that in 
case of a ferromagnetic ground state the low-temperature 
dependence of x is expected to be different. Thus for 
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TABLE II: Electron numbers n max < n < TV for which saturated ground-state ferromagnetism in the limit U — > oo exists. Note 
that for TV > 20 the size of the Hamiltonian matrix becomes very large. Therefore, for TV > 20 only a few sectors n > n max are 
accessible by numerical calculation. In particular, for TV = 24 we can examine only the sectors with up to n = 16 electrons and 
for TV = 30 with up to n = 10 electrons. Note further that for TV = 40 and TV = 50 we cannot reach the sectors with n > n max . 
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the spin- 1/2 ferromagnetic Heisenberg chain the Bethe 
ansatz yields x ^ T~ 2 (see Ref. which holds also for 
weakly frustrated chains^ 

Again we have confirmed Eq. (|6.8p by numerical cal- 
culations for finite systems. As an example some results 
for the sawtooth chain of two different lengths are pre- 
sented in Fig. H2 for t' = y/2 and t' = 1. From Fig. [HI it 
is obvious that for the sawtooth chain with a flat band, 
i.e., t' = y/2, Eq. (|6.8[) holds at low temperatures. More- 
over, the obtained Curie constants imply (S 2 )4 = 4 for 
Af = 8 and {S 2 ) A = 60 /17 fo r Af = 12 in agreement 
with calculations of Sec. IVI Al The temperature region 
where this relation is valid increases with growing U . In 
case of a dispersive lowest band, i.e., t' = 1, the quan- 
tity (S 2 ) n depends on U and n. Indeed, it has been 
observed previously^ that in the sawtooth chain with 
t' = t for quarter filling and less than quarter filling the 
Coulomb repulsion U may drive transitions from singlet 
ground states present at small values of U to ferromag- 
netic ground states present at large values of U. For the 
values of U and n considered in Fig. [12] we have singlet 
ground states, i.e., (S 2 ) n = 0, for t' = 1. Hence, contrary 
to the flat-band case we have T\(T, n, N) = as T — > 0. 
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FIG. 12: (Color online) Uniform magnetic susceptibility 
3Tx{T,n,N)/C n for the sawtooth chain (t = 1) with n = 4 
electrons and TV = 16 sites with [7 = 4 (filled triangles and 
circles) as well as TV = 24 sites with U — > oo (empty triangles 
and circles). Triangles correspond to t' = \p2 and circles cor- 
respond to t' — 1. Here we use for the normalization of the 
vertical axis C„ = 4 for Af = 8 and C n = 60/17 for Af = 12. 



VII. RELATION TO THE XXZ MODEL 

Finally we want discuss the relation between the exact 
many-electron states considered in this paper and the 
localized magnon states found for the XXZ Heisenberg 
antifcrromagnet i 26 i 29 i 30 

First we notice that the localized magnon states for 
the XXZ Heisenberg antiferromagnet on all three lat- 
tices can be also mapped onto the model of hard dimers 
on an auxiliary simple chain. However, in contradistinc- 
tion to the electron model, the localized magnons cannot 
sit in neighboring traps and there is only one possibil- 
ity to occupy a trap. Moreover, for the kagome chains 
the simultaneous occupation of a leg and a diamond or 
hexagon trapping cell by magnons is not allowed. Hence, 
we are faced with an example where the Pauli princi- 
ple leads to less constraints in comparison to hard-core 
bosonic systems. As a result, the number of states for 
the XXZ model is given by hard dimers on a chain of 
only Af sites instead of the 2Af sites as is the case for the 
Hubbard model^i^o 



Further differences can be emphasized between the lo- 
calized magnons and the localized electrons for these 
three chains. First, in the subspace with n = Af/2 
magnons the ground-state degeneracy equals 2 whereas 
in the corresponding subspace with n = Af electrons the 
hard-dimer ground-state degeneracy equals Af + 1. For 
the antiferromagnetic Heisenberg kagome chains the ex- 
tended single-magnon states (two-leg states) also appear, 
however, because of the stronger constraint only in the 
subspaccs with n = 1 and n = 2 magnons. This increases 
the total ground-state degeneracy at the saturation mag- 
netic field h = h sa ,t exactly by 2 (compare Ref. 30b). By 
contrast, for the Hubbard kagome chains the extended 
single-electron states (two-leg states) appear in all sub- 
spaces with n = 1, . . . ,J\f + 1 electrons thus noticeably 
increasing the ground-state degeneracy at jj, = 

Despite the differences between the fermionic and the 
hard-core bosonic systems stressed above, the existence 
of localized states in both systems leads to some com- 
mon features. For instance, the jump in the number of 
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electrons as a function of the chemical potential /i at 
fi = /Ug found for the electron system corresponds to a 
jump in the magnetization curve as a function of the ex- 
ternal magnetic field at the saturation field h = /i sa t seen 
in the XXZ model. Moreover, in both cases the con- 
tribution of the localized states to the partition function 
can be calculated explicitly by a transfer-matrix method, 
which leads to simple analytical expressions for the low- 
temperature behavior of various thermodynamic quanti- 
ties, such as the entropy or the specific heat. 
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VIII. CONCLUSIONS 



APPENDIX A: TRANSFER-MATRIX 
COUNTING OF LOCALIZED HARD-DIMER 
ELECTRON STATES 



In this paper, we have considered two different types 
of flat-band ferromagnets, namely, the sawtooth Hubbard 
chain (Tasaki's model) and two kagomc Hubbard chains 
(Miclke's models). For these three models we have con- 
structed the complete set of ground states for electron 
numbers n < ?i max , where n max = Af (sawtooth, odd- Af 
kagomc I) or n max = Af + 1 (even- Af kagome I, kagomc 
II) with Af = N/2 (sawtooth), Af = N/3 (kagome I), 
and Af = N/5 (kagome II). In these ground states the 
electrons are trapped on restricted areas of the full lat- 
tice. Using a transfer-matrix method, see Appendix \^ 
we have calculated the degeneracy of the ground states 
gj^(n). The ground-state degeneracy grows rapidly with 
increasing system size, and a finite residual entropy per 
site remains in the thermodynamic limit. Moreover, we 
have calculated exactly the contribution of the highly 
degenerate ground-state manifold to the partition func- 
tion. The low-temperature thermodynamics around a 
particular value of the chemical potential /io (/io = 2t for 
all three models) is dominated by these trapped ground 
states leading to a low-energy scale separated from the 
usual energy scale determined by the band width. In the 
thermodynamic limit N — > oo all models exhibit identical 
thermodynamic behavior in this regime which is analo- 
gous to that of classical one-dimensional hard dimers. 

The trapped ground states lead also to particular mag- 
netic behavior including ground-state ferromagnctism as 
well as paramagnetic behavior. 

Moreover, with this study we have (i) illuminated rela- 
tions between Tasaki's and Mielke's flat-band ferromag- 
nets, and (ii) the relations between frustrated quantum 
Heisenberg antiferromagnets and the Hubbard flat-band 
ferromagnets. In spite of some similarities in the math- 
ematical description of both correlated quantum lattice 
systems owing to localized one-particle states, the elab- 
oration of a comprehensive theory for the Hubbard flat- 
band fcrromagnct in higher dimensions using some ideas 
from localized-spin systems remains an unsolved problem 
and calls for further efforts. 



The ground states of electrons localized on trapping 
cells (V- valleys, diamonds or hexagons, respectively) for 
U > can be counted using a transfer-matrix method^ 
In Ref. Qj] we have associated two sites to each trapping 
cell and used a hard-dimer mapping for this counting. 
Here we will present a direct solution of the same problem 
using a 3 x 3 transfer matrix. 

Recall that the ground states of the considered sys- 
tems, i.e., the sawtooth chain or the odd- Af kagome 
chain I with /i around can be obtained by populat- 
ing the trapping cells with electrons with spin-up and 
spin-down according to the following rules: 1) Each trap 
may be empty, occupied by one spin-up electron or oc- 
cupied by one spin-down electron. That is, for each trap 
j = 0, 1, . . . ,Af— 1 we have three trap states, Sj = 0, f , J,. 
A ground state of the chain can be thought of as a certain 
sequence of the trap states. 2) Moving along the succes- 
sive traps (the choice of the first trap is totally arbitrary, 
e.g., we may take j = 0) we must allow only one sequence 
of two trap states corresponding to the neighboring cells 
being occupied by differently polarized electrons. This is 
only a convention for the correct counting of the number 
of the ground states. (Note, however, that the convention 
does not work for the number of electrons n = Af yield- 
ing only 2 states instead of correct number Af + 1.) For 
instance, let us allow the sequence of trap states Sj =|, 
sj + i =1 and forbid the sequence of trap states Sj =1, 
Sj+l =t 

These rules can be encoded with a transfer matrij*^ 
T(0,0) T(0,T) T(0,i) 

T(t,o) r(u) r(u) 

T(|,0) T(U) T(u) 



ill 

z z z 
z z 



(Al) 



with the activity z = exp(— e\/T). Now we can 
write down the contribution of all allowed sequences 
So, . . . , Sjv"-i to the grand-canonical partition function 
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Strap (Zj -AO <|5.2p as follows 



Strap [Z, Af) 



TrT^ 



so 



(A2) 



The hard-dimer computation yields the alternative rep 
resentation Strap (<£) -A/") 
•^/l /4 + exp x are the eigenvalues of a 2 



A 2AA + A 2AT where a ± = 1/2 ± 



2 transfer 

matrix.— Since £± = A± and £o = 0, the two expres- 
sions are in fact equivalent. 

We turn now to the canonical description. More specif- 
ically, we consider n < Af trapped states on a (periodic) 
chain with Af cells. We are interested in the canonical 



partition function of such a system Z(n,Af) which counts 
the number of spatial configurations of n trapped elec- 
trons. Using the relation between the canonical partition 
function and the grand-canonical partition function 



(A3) 



n=0 



we immediately find that 

i crs trap (z,A0 



Z(n,Af) 



dz r - 



(A4) 



z=0 



After simple (but becoming tedious as n increases) cal- 
culations we get 



Z(l,Af) = 2Af, 
Z(2,Af) = M(2JV-3), 

Z(3,Af) = ^2AA(4AA 2 - 18AA+20), 

Z(4,Af) = ^2Af(8Af 3 -72Af 2 + 2UAf -210), 

Z(5,Af) = ^-2AA(16AA 4 - 240A^ 3 + 1340AA 2 - 3300AA+3024), (A5) 
5! 

Z(6,Af) = ^-2M(32M 5 - 720AA 4 + 6440A/" 3 - 28620A/" 2 + 63188A/" - 55440), 
6! 

2(7, AT) = ^2Af(64Af 6 - 2016A/" 5 + 26320A/" 4 - 182280A/" 3 + 706216A/" 2 - 1451184A/" + 1235520), 

2(8, AO = ^2A"(128AA 7 - 5376A/" 6 + 96320A/" 5 - 954240A" 4 + 5645192A/" 3 - 19941264A/" 2 + 38943000A/"- 32432400). 
8! 

The data shown in Fig. [3] are based on Eq. (|A5|) . 



Let us calculate the (not normalized) grand-canonical 
correlation function (S§Sj) z , where 5| is the z- 
component spin operator of the trap j and the subscript 
z denotes the activity. Defining a matrix 

/ \ 
S = \ (A6) 

\o 6 -\) 

one passes again from the sum over hard-dimer ground 
states for a fixed number of electrons n and the sum over 
n < Af to the sum over sq = 0,1,1, s i = OjTjIj ■ ■ ■ ■> 
sjv-i = 0, t, 4- As a result 

(SJS/), = Tr (STJST^') . (A7) 

Using a MAPLE code we can easily compute (S§Sj) z 
according to Eq. (|A7j) for sufficiently large systems (up 
to Af = 256); the resulting expression for {S^Sj) z is a 
polynomial with the powers of z from 2 to Af. 



Turning to the canonical description we use the rela- 
tion 

N 

(S%S?) Z = z n Z(n,Af)(S$S?) n (A8) 

n=2 

[obviously (<5o-S|) n =o = (S^Spn^i — do not enter the 
right-hand side of Eq. (|A8|) ] to derive 



(A9) 

2=0 



1 ePiSZS* 
Z(n,Af)(S* S?) n = - 



77 



! dz r 



Thus the coefficients associated with the corresponding 
powers of the activity z in the right-hand side in Eq. (|A7[) 
yield the required quantities (S^Spn. 

Finally, we mention that the same results for (S 2 ) n 
can be obtained within the hard-dimer picture. Since 
each hard-dimer ground state can be represented in terms 
of the occupation numbers of hard dimers on a sim- 
ple chain, i.e., it is enumerated by a set of hard-dimer 
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occupation numbers no, ni, . . . , n^U-ii n j — 0,1, and 
S z = (1/2) E^T^-iy'nj, we have 

(S z \ = y E (-l) 9 (^on 9 )„, (A10) 

9=0 

where (. . .}„ in the right-hand side in Eq. (|A10|) stands 
for the (normalized) average over spatial configurations 
of n hard dimers on a simple chain of 27V sites. To 
find a density-density correlation function at distance 
q for one-dimensional hard dimers in the canonical en- 
semble (n, 27V), it is convenient to calculate first a 
density-density correlation function at distance q for one- 
dimensional hard dimers in the grand-canonical ensemble 
(z,27V), 

(n n q ) z = Tr (ND«ND 2A ^) , 

■>-Utf). »-(!:)■ <*»> 

which gives the required (non q ) n after inverting the re- 
lation (n n q ) z = J2n=2 z n Z(n,N)(n n q ) n . 



APPENDIX B: LINEAR INDEPENDENCE OF 
TRAPPED ELECTRON STATES 

We wish to clarify whether the set of ground states 
constructed in Sees. IfV Al and IIVBI for given n < n max 
is linearly independent. An affirmative answer for hard- 
dimer states comes from Ref. H^. For the sawtooth chain 
and the kagome chains the localized n-clcctron states 
(n = 1, . . . ,7V) are linearly independent, which is con- 
nected with the fact that for all three lattices there are 
sites which are unique to each cell (isolated class in the 
nomenclature of Ref. [48h . 

We can use the same arguments for the set of single- 
electron states which consists of localized states on a di- 
amond/hexagon plus one two-leg state. For this purpose 
the two-leg state (|3.5[) / (|3.9p in the set of states can be 
replaced by the state localized along one (e.g., lower) leg 
only. As a result we again are faced with the case when 
there are isolated sites [the sites belonging to another 
(upper) leg] that yields linear independence of the con- 
sidered set of single-electron states and thus of n-electron 
states (n= 1, . .. ,7V + 1)M 



The Hubbard Model — A Reprint Volume, edited by 
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